Universal width distributions in non-Markovian Gaussian processes 
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We study the influence of boundary conditions on self-affine random functions u(t) in the interval 
t/L G [0, 1], with independent Gaussian Fourier modes of variance ~ l/q a - We consider the prob- 
ability distribution of the mean square width of u(t) taken over the whole interval or in a window 
t/L G [x,x + 8]. Its characteristic function can be expressed in terms of the spectrum of an infinite 
matrix. This distribution strongly depends on the boundary conditions of u(t) for finite 5, but we 
show that it is universal (independent of boundary conditions) in the small-window limit (6 — > 0, 
5 <C min[x, 1 — x]). We compute it directly for all values of a, using, for a < 3, an asymptotic 
expansion formula that we derive. For a > 3, the limiting width distribution is independent of a. It 
corresponds to an infinite matrix with a single non-zero eigenvalue. We give the exact expression for 
the width distribution in this case. Our analysis facilitates the estimation of the roughness exponent 
from experimental data, in cases where the standard extrapolation method cannot be used. 



I. INTRODUCTION 

In nature, random processes model interfaces and surfaces [1] [2], turbulent flows [3], erratic time series[4] and many 
other systems. In the most simple setting, these random processes correspond to a function of a scalar variable u(t), and 
are characterized by a probability distribution V[u(t)] that stems from an equilibrium problem or a non-equilibrium 
process. The probability distribution V[u(t)] is often unknown, and generally inaccessible to exact analysis. In many 
cases, approximate methods must be brought to bear on these problems. 

One of the most successful approaches in the field of random processes is the Gaussian approximation. It consists, 
schematically, in developing u(t) on a general Fourier basis, 

00 

u(t) ~ ^ a n exp (ig n i) , (1) 

n— — oo 

and in assuming that the coefficients a n = a*_ n are independent Gaussian random variables. The probability distribu- 
tion V[u(t)] = V[{ai,a,2, ■ ■ ■}} then factorizes into a product V ~ p Gauss — V {a\)V {a-i) . . . , where V(a n ) is a Gaussian 
with zero mean and variance <j\. The idea behind this representation is to approximate the translation- invariant action 
of a complicated interacting problem by a quadratic (Gaussian) action which yields the same two-point correlation 
functions as in the original problem. In the particularly important case of self-affine (critical) systems, the only length 
scale present is the system size. The scaling of the variances with the modes q n can then be described by a single pa- 
rameter a with a\ cx l/n a . In many non-trivial problems, the Gaussian approximation is in outstanding quantitative 
agreement with the full theory [5, 6], and reproduces very well even the higher-order correlation functions. 
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FIG. 1: Gaussian functions u(t) in the interval t/L 
C = 0.75 (a = 2.5) (full periodic series). 



£ [0,1], corresponding to a probability distribution - P Gauss [u(t)] with 



Because of its decoupled Fourier modes, the Gaussian approximation is considerably simpler than the full theory 
and the function u(t) can easily be generated through eqn (1) from Gaussian random numbers {a n } (see Fig. 1). 
However, the price to pay for this simplicity in Fourier space is to have non-trivial long-range correlations in real 
space. As we will discuss in detail, the real-space action contains generalized derivatives which, in the sense of the 
Ricmann-Liouville derivative, can be expressed as an integral convolution with a long-range kernel [7]. For non- 
integer values of a/2, the real-space Gaussian action is non-local and the geometrical properties of the function u(t) 
are intricate. 

For a = 2, the Gaussian approximation corresponds to the notorious random walk. In this case the real-space action 
is local. It defines a Markovian evolution (the value of u at t + dt depends only on the one at t). For this reason it is 
possible to determine in detail its geometric properties (for a recent review sec [8]). For a ^ 2, instead, the process 
u(t) is non-Mar kovian: In the case a — 4 (the driven curvature model [1, 9]), the evolution of the derivative u(t) 
is Markovian, but not the one of u(t) itself. For not-integer a/2, the non-Markovian properties reflect the non-local 
character of the action. This implies that, generally, memory effects influence the shape of u(t) and calculations, even 
within the Gaussian approximation, arc difficult. In particular, the persistence exponents and the distribution of the 
extreme remain unknown [10, 11]. On the other hand, the boundary conditions influence the statistical properties of 
u(t) in a non-trivial way for all t G [0, L\. Understanding these effects is important because the periodic boundary 
conditions for u(t), which are commonplace in theoretical calculations and in numerical simulations, are usually not 
realized in experiments [12, 13]. 

As an example of a fundamental geometrical quantity sensitive to the boundary, we consider the mean square width 
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which is relevant both from the theoretical and the experimental points of view. For a self-affine random process, the 
mean square width is itself described by a non-trivial probability distribution P(ui 2 ), which is neither a delta- function 
nor a Gaussian. 

In this work, we study the influence of boundary conditions on the width distribution of V uss [w(i)] for self-affine 
functions characterized by variances scaling with the single parameter a. In fact, the schematic Fourier representation 
of eqn (1) can be rendered explicit in a number of ways in order to accommodate the boundary conditions. First it 
can be expanded into a full periodic series, 



u(t) 



OO 



/27m 
, a n cos I —j—t 

n=l ^ L 



(2im 
— ^— i ) (full periodic series), 



(2) 



where a n and b n are independent Gaussian random numbers of variance l/n a . This implies that the function, 
which has zero average value, and all the derivatives are periodic, if they exist (u^(t = 0) = u^ k \t = L) for 
k = 0, 1, • • • < a/2,). It is also possible to expand the Gaussian function u(t) into a sine Fourier series, 



u(t) — ^ s n sin (J-jrtJ (sine series), 



(3) 



n=l 



again supposing that the s n are independent Gaussians with variance ocl/n a . In that case, the function u(t) vanishes 
at t = and t = L. By a uniform shift of the function, it can be made to have zero average value, as for the full 
periodic series. However, all the existing even derivatives vanish (u^ k \0) — u^(L) = for k = 0, 2, . . . < a/2) for 
the sine Fourier series. The full periodic series and the sine series are equivalent for the random walk (a = 2) [14, 15], 
but they differ for all other values of a. 

Finally, the function u(t) can also be expanded into a cosine Fourier series, 



oc 

u(t) — ^ c„ cos (-jrtj (cosine series), 



(4) 



n=l 



3 



again with a 2 {c n ) oc l/n a . In this case, the random function is not forced to satisfy u(0) = u(L), and for this reason 
it has been used to study free random walks [8]. Analogously to the sine Fourier series, all the odd derivatives of u(t) 
must vanish at the boundaries, if they exist (u^ k \0) = vS k \L) = for k = 1,3, ... < a/2). This paper shows that 
these additional constraints strongly influence the geometry of the function u(t) for a^2. 

The above three series representations correspond to the same real-space action <S[u(t)] and its associated Gaussian 
probability distribution -pGa USS 
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S[u(t)\ = -] 



dt 



d a ' 2 u{t) 

Q t a/2 



pOauss [ u ( t )] =exp {_5[ u ( t )]}. 



(5) 



/2] e iqt = 



For non-integer values of a/2, the generalized derivative in eqn (5) is defined in momentum space: [d a ^ 2 /dt a 
{\q) a / 2 e lqt (the derivative can also be defined in real space [7]). Using this definition we can compute the variances of 
the Fourier coefficients: 



7> Gauss [K,6„}] = Hrr^p 



" 1 2 a - 1 iT a n a 
For the sine and the cosine series, the variances are 
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(full periodic series). 
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(6) 



(7) 



This choice leads to analogous expressions for the probability distributions -p Gauss [{s n }] and •p Gauss [{c„}], respec- 
tively. 

The variances in eqns (6) and (7) scale with the system size as oc L a ~ 1 = L 2 ^ where £ is the roughness exponent. 
This exponent characterizes the main geometric properties of a self-affine system. 
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FIG. 2: Gaussian functions u(t/L) corresponding to the full periodic and the cosine series. Window boundary conditions 
correspond to picking out a piece of u(t) in the interval t/L € [x, x + S] and shifting u(t/L) such that the window mean value 
vanishes. 

In this paper, we first consider Gaussian functions on the entire interval [0,L]. We compute the average of the 
width distribution for the three series and for general values of a (Section II A). The average of the width distribution 
strongly depends on the boundary conditions (for a / 2). These results are then generalized (Section II B) to the 
case when the function u(t) is restricted to a window of width S, in the interval t/L e [x, x + 5]: 



w 2 {x,5) = - 



dt u 2 (tL) 




dt u(tL) 



(8) 



These window boundary conditions are closer to the experimental situation than those realized by either the full 
periodic or the cosine series (see Fig. 2). We write (Section III) the characteristic function of the width distribution 
in terms of the eigenvalues of an infinite matrix which depends on the basis functions of the Fourier series, and on the 
window parameters {<5, x}. Practically, to compute the width distribution for any of the series, in an arbitrary finite 
window, it suffices to compute this characteristic function from the eigenvalues of a finite-rank approximation of the 
above matrix, and to perform an inverse Fourier transform. It is easy to see that the width distribution depends on 
the size of the window and on the choice of boundary conditions. In the small-window limit (5^0 and far from 
the boundaries S <C min[x, 1 — x]), finite-rank approximations to the above matrix cease to be accurate (the limit of 
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FIG. 3: Mean square width as a function of window width, for different values of a, for the full periodic series (left) and the 
cosine series (right, with x = |(1 — 5)) (from eqns (11) and (16)). 

rank N — > oo docs not commute with the limit 5 — > 0, for a < 3), but we are able to write all the cumulants of the 
width distribution distribution directly in this limit (Section III C) via a subtle asymptotic expansion in powers of 5. 
We prove that in the small-window limit the width distribution becomes independent of the boundary condition. We 
also show how to practically compute the (a-dependent) universal width distribution directly in this limit. For a > 3 
(Section IIIC2), the problem simplifies. The characteristic function in the small- window limit then corresponds to a 
matrix with only one non-zero eigenvalue, and the corresponding width distribution no longer depends on a. We give 
its explicit form. Finally, we compute the logarithmic corrections to these asymptotic results for odd-integer values 
of a. In Appendix A and Appendix B, we provide technical details of our calculation. 

In a previous paper on the same subject [16], we already studied the second moment of the width distribution, 
and presented arguments for the universality in the small- window limit. The more complete, and more concrete, 
calculations of the present paper rely on the representation of the characteristic function in terms of the spectrum of 
a matrix, which was not contained in [16]. 

II. AVERAGE WIDTH 

In a self-affine (critical) system, the length L of the total interval is, as mentioned, the only characteristic length 
and the average value of the total width scales as (-u^) oc ((. . .) denotes the ensemble average). 



A. Average width for 8 = 1 



Fixing <5 = 1 and x = in eqn (8) and integrating over t yields 
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(full periodic series) 
(sine series) 
(cosine series). 



From the Gaussian probability distributions given in eqn (5)-eqn (7), the ensemble averages are 

(full periodic series) 
(27r)° 1 ^( a + 2) (sine series) 



(w 2 ) 
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£C(a) - 2 2 " +2 - 



l£C(«) 



(cosine series), 



where C( x ) — E^Li ^/ nX is the Riemann Zeta function. For all a, the average width of the full periodic series is 
smaller than that of the cosine series, as is quite natural. For a = 2, the average widths of the full periodic series 
and the sine series coincide, because the two differ only through boundary conditions for the derivative of u, that a 
Markovian process is insensitive to. 
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B. Average width for S < 1 



We now compute the average width for S < 1 for the full periodic series, the cosine, and the sine series. The width 
u, 2° r (^) i s independent of the origin x but the x-dependencc of the average width cannot be neglected in the other 
two cases. For the full periodic series, one obtains from eqns (2) and (8) 

oo 

wf 1 ^) = E a n a m C nm (S) + &nb m Inm{3) ~\~ b n b m S nm 

(5), (9) 

n.m— 1 

where the coefficients C nm (S) are given by elementary integrals: 

px+6 -i rx+5 rx+8 



Cmn = t / di cos(27rmt) cos(27T77i) — — / At cos(27rmt) / At cos(27mt). 
5 Jx ° 2 J x J x 



(10) 



Analogously, S nm and I nm can be expressed in terms of sine-sine and cosine-sine integrals. For the full periodic 
series, these coefficients are naturally independent of x. 

Equation (9) allows to compute for one given sample. Integrating over the Gaussian Fourier components 
{a n ,b n }, we get 

/ per/c-w L a 1 sr^ C nn + S nn 1 — cos(27rn<5) . . 

K («)> = E na ■ where C ™ + S ™ = 1 - 2 (nnSV ■ 

n— 1 v y 

The sum in eqn (11) is easily evaluated for finite S (see Fig. 3). The limit of this mean value for 5^0 cannot be 
obtained by a naive Taylor expansions of each term in this infinite sum in eqn (11), because it is not uniformally 
convergent in the interval S <G [0,1] for a < 3 (in the limit 5 — > 0, the terms of the sum behave as (C nn (S) + 
S n n{5))n~ a ~ <5 2 n 2 ~ Q , producing a diverging series for a < 3). As we will discuss in detail later, the higher cumulants 
of the width distribution are given by multiple infinite sums which present the same pathology as the sum in eqn (11). 
To sum the series in the limit 6 — ► 0, we have derived a very useful expansion formula which has the same structure 
as the Euler-Maclaurin formula: 

v« = ^ r dt L v a^) + f r/ M (0 )M d2) 

^ n a J \ t a ^ ml ^ ml V ' 

n=l u \ m=0 / m=0 

where \a\ is the integer part of a. Equation (12) holds inside the convergence radius 5 — 1 for all the quantities 
considered in this work (the formula is proved in Appendix A). For analytic functions f(z) and non-integer a, the 
first term on the right can be interpreted as the naive limit of the sum as S — > 0, with t = nS, whereas the second 
term contains the Taylor expansion of f(n5) around zero. For integer a, the singularity of £(1) generates additional 
logarithms (see Appendix A, again). Using the expansion formula of eqn (12), one arrives at: 

« er («5)) 2— 1 a a + 2) sa _ 1: 4 ~ 1 , n+1 C(«-2n) 



L"- 1 ((-a -I) 7r" +2 ' ( 27r ) a ~[ (2n + 2)! 

To take into account the logarithmic corrections for odd integer a we must use eqn (A6) instead of eqn (12). The 
final result writes: 

{wpW = ^ M« + 2) ~ yrJ) x 4 V (-1)" +1 #^S(2^) 2 " (« odd integer), 
L"- 1 v 7 27r(a + l)! (27r)« v ; (2n + 2)P 7 v 6 " 

where ^o(-z) is the digamma function. An expansion analogous to eqn (13) appears in the correlation function 
governing the density of zero-crossings of a Gaussian function [10]. For integer even a, the series in eqn (13) is finite 
(because ((x) = for x — {—2,-4, ...}). The mean square widths for the periodic random walk (a = 2) and the 
driven curvature model (a = 4) then have compact expressions: 



«"•(£)) _ fjjS-^S 2 fora = 2 

L"-l -\ T 1 I< 52__1_ 5 3 + _1_ (5 4 fora = 4 
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For the cosine series, we obtain, from eqns (4) and (8), the mean squared width 

oo 

w c 2 os {x,5) = c n c m D nm (x,S), (14) 



n,m=l 

where the coefficients are given by a symmetric matrix 



rx+5 y rx+8 rx+8 

D mn (x,S) = - / dt cos(irmt) cos(irnt) — — / dt cos(7rmt) / dt cos(7rni). (15) 

° Jx ° Jx Jx 

D nm is the overlap matrix of the basis functions in eqn (4) on the interval [x, x + 5]. For 5 = 1, the basis functions 
are, by construction, an orthonormal set, and D nm = \5 nm . 

Integrating over the Gaussian variable {c„}, the average width (u>2° s ) becomes 

(^,5)) = 2 ^±^^1. (16) 



n=l 



The sum in eqn (16) is again easily evaluated, for any value of x and 5. The behavior of (w2° s (x, 5)) as a function of 
5 is shown in Fig. 3 for x = |(1 — 5). For small 5 we notice that (w2° s {x, 5)) oc <5 Q_1 for a < 3 and (w2° s (x, 5)) oc 5 2 
for a > 3. A special case is again a = 2 where we have (w2° s (x, 5)) oc 5 for all S. 



Small-window limit 



In the following, we extract the universal behavior of {w2° s (x, 5)} for a small window [a;, x+S] far from the boundaries 
(for 5 <C min(a;, 1 — x)). For a < 3, eqn (16) can be expanded as 
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where the function D(t, x/5) is defined through D(nnS, x) — D nn (x). We write D(t, x/5) as a sum of two terms: 

D(t, x/5) = h(t, x/5) + 1 [C(t) + S(t)} . 

The function h(t,x/5) contains all the x-dependence of D(t,x/S): 

sin(2i) 



h(t, - 5 ) = cos(2i(| + l -)) 



At 



^sm(-) 



For i/i5> 1, h(t, x/5) oscillates rapidly and gives a vanishing contribution to the integral in eqn (17). In particular, 
we can show that 

(w?*(x,5)) 2— 1 aa + Vsc-i 



lim 

x/S^QO 



r a-l 



C(-a-l) tt«+ 2 



■0[exp(-x/6)] . 



Comparing the above expression with eqn (13), we conclude that, for a < 3, the dominant contributions to (w2° s (x, 5)} 
and to (w2 Cr {5)) coincide for small 5 (compare also with Fig. 3). This conclusion can be extended to the sine series, 
where the explicit x-dependence is again due to oscillatory terms that vanish in the small- window limit. 

As mentioned before, the series in eqn (16) is uniformly convergent for a > 3, so that one can take the limit 5^0 
for each mode and then do the sum. Using the expansion 



D nm {x) = -^nmir 2 sin^mx) sm(Trnx)5 2 + 0(5 4 ) 



in eqn (16), we obtain 



(w C 2 OS (x,5)) 



r a-l 



y~v sin (nirx) 

.71=1 
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a-2 



"5 2 + 0(S a -\5 i ). 



(18) 



(19) 



For x = and x = 1 the 5 2 term vanishes, while for < x < 1 it is smaller than the corresponding term for the 
periodic case. This behavior is consistent with the fact that, for a > 2, the cosine series imposes vanishing derivatives 
at the end points. Hence these boundary conditions force (w™ 8 ) to be smaller than {w^} for all < x < 1. The sine 
series gives a result analogous to the one for the cosine series. It suffices to replace the sine in the sum of eqn (19) by 
a cosine. 
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III. WIDTH DISTRIBUTION 



The scaling behavior of the average width gives access to the value of a and thus to the roughness exponent 
Q = (a — l)/2. However, the value of this exponent only depends on the 2-point correlation functions, and captures 
no finer geometric properties of the function u(t). To discriminate between a Gaussian and a non-Gaussian function, 
one must have access to higher cumulants, as they are contained in the sample-to-sample fluctuations of the two-point 
correlation functions. The distribution V(w2) has been used to analyze numerical and experimental data [13, 20- 
23]. The width distribution also allows to estimate the roughness exponent from experimental data which are not 
sufficiently good to plot (u>2(S)) vs S over several orders of magnitude in 5. 

For general values of a, the width distribution has been computed for the entire interval (6 = 1, x = 0), for the full 
periodic series, [17-19], and for the cosine series [13]. The underlying simplification with respect to the calculations in 
the present paper is that the matrix of coefficients for 8 = 1 satisfies I nm = while C nm , S nm and D nm are diagonal, 
as evident in eqns (10) and (15). In order to study window boundary conditions (but also for the sine series), the 
previous framework must be generalized to non-diagonal matrices. 

The width distribution can be obtained from the symmetric matrices 

A™ = 2a c n ° s D nm a™. (20) 

and 

/ 9rrP cv d rrP° r rrP or T irP or \ 

I rrP CI T rrP QI 9rrP or <? iTP or I ' ^ ' 

\ u n 1 mnO m n Onm<-> m J 

respectively. Concretely (see eqns (9) and (14)), individual realizations of the mean square width distribution are 
generated by multiplying the matrices A in eqns (20) and (21) by vectors of normal distributed Gaussian variables. 
This allows to obtain V(z = u> 2 / (w^)) approximately through direct simulation (see [15]). 

On the other hand, one can compute all the moments of the distribution from contractions of a given matrix A 
(which can stand for A cos , A sm , or A pcr ). For example, the second moment of the width distribution is given by 



= <»2 2 ) = j£ A™ A™ C) - ^Tr L4 c °f + 



-TtA cos 
2 



(22) 



Similar expressions exist for higher moments of the width distribution. As shown in Appendix B, the cumulant ki of 
the rescaled width distribution can be expressed in a simpler way than the moments, as a trace of the matrix A taken 
to the Zth power: 

2 {102)' 1 1 

This is already apparent in eqn (22), where fi 2 = { w 2} 2 («2 + «a). The cumulants are thus given in terms of the 
eigenvalues {Ai, A2, . . . } of the matrix A as TrA 1 = J^k 

(*-!)! EfcAj 

Kl= 2 (j: k x k y (23) 

The cumulants in eqn (23) yield an explicit formula for the cumulant-generating function ^(s/ (^2)), where 

00 00 1 

^) = EE^ = E^ 1 °g( 1 -^)- 

fe=l 1=1 k 

The characteristic function f(s) = exp(^(is)), the exponential of the cumulant-generating function Vf'(is), is given by 

/( s ) = iT^r^=- ( 24 ) 

k V 1 - lAfeS 

The spectrum of the matrices under consideration is such that the infinite product in eqn (24) is uniformally convergent. 
The associated width distribution is recovered through an inverse Fourier transform, 

V{z = w 2 / (w 2 )) = -L / ds cxp(-izs)/( S ), (25) 

which can be obtained by straightforward Riemann integration because the branch points s k = i/A fe are away from 
the real axis. 
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mean square width (rescaled) z mean square width (rescaled) z 



FIG. 4: Rescaled probability distributions for the mean square width for £ = 0.75 (a = 2.5) for the cosine series (left) and the 
full periodic series (right), obtained from eqns (24) and (25). The result of direct simulations for 8 = 1/16 (from eqn (14)), and 
the solution in the limit 5—^0 (from the matrix in eqn (32)) are also shown. Sizes of matrices are N = 512 and N = 1024. 



A. Width distribution for the entire interval (8 = 1) 



As mentioned above, in the case 5=1, the matrices 

l2 



are diagonal and the computation of the cumulants from eqn (23) is direct: 



Kl = 



(I — l)\((la) I '((a) 1 (full periodic series) 
(21 - 2)!!C(/a)/C(a)' (cosine series). 



The cumulants for the sine series are more complicated because the corresponding matrix ^4 sm presents non-diagonal 
terms coming from the non- vanishing mean dt sin (nnt/L) for odd values of n. The second cumulant, for 

example, is given by 

sin _ 2C(2a) + 2- 2q+1 (16tt- 4 (2 q + 2 - l) 2 ((a + 2) 2 - 8^ 2 (2 2tt + 2 - l)C(2a + 2)) 
K2 ~ [((a) - 2-«+ 1 7r- 2 (2«+ 2 - l)((a + 2)] 2 ' 

K2 n agrees with only for a — 2, as expected. In the case of the full periodic and the cosine series, the characteristic 
functions assume a simple form 



/(*) 



n fe (l - iL a - 1 (2Tr) a k- a s)- 1 (full periodic series) 



Ik 

11^(1 _ iL a ~ 1 n a k~ a s)~'5 (cosine series). 



As already discussed in [19], the two-fold degeneracy of the spectrum of A per for 5 = 1, evident in eqn (26), yields a 
characteristic function f(s) with simple poles on the imaginary axis. This simplifies the inverse Fourier transformation. 



B. Width distribution for windows (finite <5) 

The expressions allowing to recover the width distribution from the eigenvalues of the matrices A cos ,^4 per , and 
A sln remains valid for intervals S < 1, even though the computer must now be used for finite approximations of 
these infinite matrices, which are no longer diagonal. These calculations can be easily checked by direct simulation, 
as mentioned before. The outcome of this analysis is shown in Fig. 4 for the case a = 5/2. It is evident that the 
width distribution changes with the sample size for the value of a chosen. Furthermore, the direct evaluation of the 
rescaled width distribution for finite but small 5 suggests that this distribution becomes universal (independent of the 
boundary conditions) in the limit 5^0. This scenario is confirmed by comparing the evolution of the (normalized) 
spectra of A per and A cos , as shown in Fig. 5. 
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FIG. 5: Eigenvalues {Ai, A2, . . . } for the cosine and the full periodic series as a function of window size 5 ( a — 2.5 (£ = 0.75), 
for the cosine series, x = |(1 — 5)). The spectrum in the small-window limit 5^0 (crosses), is obtained in Section IIIC. 
Eigenvalues are normalized as J^ fe A fe = 2. 



We then compute V(z = w 2 / (w 2 )) directly in the limit 8 — > 0. However, for small 5, increasingly larger matrices 
A pcr and A cos must be considered because of the non- uniform convergence of the traces of these matrices for 8 € [0, 1]. 

In Section IIIC we determine directly the width distribution at the limit 8 — > 0. In particular, analogously to the 
computation of asymptotics of the average width (Section II C), the expansion formula eqn (12) is proven useful to 
compute the cumulants of the distribution in the small-window limit for a < 3 (Section IIIC 1). This analysis serves 
two purposes: it proves the universality of the width distribution and provides a high-precision method to compute 
it directly in the limit 8 — > 0. 



C. Width distribution in the small-window limit 



In Section II C we computed the mean square widths (w 2 el (8)} and (w2° s (x, 8)} in the small-window limit. We now 
determine the rescaled width distribution V(z) in this limit. 



1. Small-window limit (a < 3) 



The computation of ki as a (^-expansion is along the lines of the previous determination of (102). For a < 3, we 
handle the non-uniform convergence of the traces in eqn (23) for 8 — > using the expansion formula eqn (12). The 
leading contribution to m, I = 1, 2, . . . , is given by multiple integrals: 

l(a-l)(] _ I \| roo poo l ji 

«f" = A -/ •••/ Y[-^Tr[A^(t 1 ,t 2 )A^(t 2 ,t 3 )---A^(t l ,t 1 )] + 0(8 3 - a ) (27) 

^w 2 j a j n=1 t n 



nr(x/8) 



2wi 



ndt n 
fa 



fa 

n=1 K 



A™{H,t 2 , - 5 )A cos (t 2 , t 3) - 5 ) ■ ■ ■ A cos (ti,ti, |)] + 0(8 3 - a ). (28) 



In the above equations, w 2 is the prefactor of the leading term in the expansion of (w 2 (S,x)), given in eqn (13): 

2-"- 1 C(a + 2) 



w 2 = 



C(-a-l) 7T«+ 2 



which is independent of the boundary conditions. A cos (t n , t n+ \,x/8) and A per (t n , t n+ \) are obtained by replacing nnS 
with t n in A 0053 and in A pel . We now verify for each cumulant 



4 os (x/«5 -» 00) -» K f cr + 0( e -^). 
Extracting the a;/(5 independent part of «^ os one obtains 



' 2^2 



TT^Tr 

11 +a 



A cos (t 1 ,i 2 )i cos (t2,i3)---i cos (ii^i) 



+ oscill. terms, 



(29) 



(30) 
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FIG. 6: Cumulant H2 of the rescaled width distribution for a = 2.5 from a naive discretization of the integral in eqn (30), and 
after the change of variables in eqn (32) (Riemann integration in the interval t G [0, 15], using A = U+i — U). 



where A s is now independent of x and given by 

4£,«-i 



A cos (t,t') 



2(t-t') hm 2 



\ 2(t+t') Sm 2 

Now it is straightforward to verify that the matrices A pcT and A cos satisfy 



t+t' 



W S ^ n 2 Sm 2 2(t+t') 
1 • t • t' 1 
it 7 Sm 2 Sm ~2 2(t-t') 



sm 
sin 



t+t' 

2 

t-t' 
2 



tF sin I sin T 



(31) 



Tr [i cos (ti, i 2 )i cos (t 2 , is) • • • i cos (ii, ti)J = Tr [A^(h/2, t 2 /2)A^(t 2 /2, t 3 /2) • • • ^ per (^/2, . 

This establishes the validity of eqn (29). 

Closed analytic expression have not been obtained for the above integrals. Actually the problem of computing the 
integrals in eqn (30) reduces to the solution of a homogeneous Fredholm equation of the second kind: 



^t,t') ij9j (t') = Affi(t), 



where the 2x2 matrix A cos (t,t') is a compact kernel with a discrete set of eigenvalues Xk converging to 0, which 
encode all the information on the cumulants Ki, in the same way as for the discrete case (eqn (23)). 

The spectrum of the kernel A cos (t 7 t') is most easily obtained by discretizing the variables t and t' on an equally 
spaced grid with tk = A(fc — \), k = 1, . . . ,N, with an upper cut-off for the integrations. However, the singularities 
of the integrands makes the convergence rather slow (~ t^ a for t n — > (with all other variables kept finite)). The 
divergence in the integrals at small t is eliminated by a standard change of variables: an integral of a function diverging 
as 1/i 7 for t — > can be written as 



[ 



dt f(t) 



1 



1-7 



dt t^/^-^f 



1/(1-7) 



The integrand on the right is constant for small t. Concretely, the change of variables in A cos leads to a matrix 

1 



A cos (t,t') = t5b^os( f l/(3-a) )t /l/(3-a)^/£r^ 



a 



(32) 



which can again be discretized. The characteristic function of the width distribution is computed from the spectrum 

of the 2N x 2N matrix A cos as discussed before. To show what is gained by rescaling the matrix A in eqn (31), 
we have computed the second cumulant k 2 on an equally spaced grid with both versions (see Fig. 6). The rescaled 
matrices converge exceptionally well with the discretization parameter A. In Fig. 5, the spectrum of the rescaled 
kernel (for 5 — > 0) is compared to the spectra of ^4 pcr and A cos for finite 8. Only in the special case a = 2 is the 
spectrum of the cosine series independent of 5, because of the Markov-chain property of the random walk. For all 
a^2, the spectrum, thus the cumulants of the width distribution, depend of 6, but satisfy 

«T(^ °) = *i os (^ °) * K i os (s = i), 

in a way illustrated in Fig. 5. As was already shown in Fig. 4 for a = 2.5, the convergence of the width distribution 
for finite 6 towards the asymptotic width distribution is quite fast, even though the rate of convergence depends on 
the value of a. 
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(cosine) S (full period.) 

FIG. 7: Eigenvalue spectrum of ^4 cos and A pcT for a — 3.5 (( = 1.25). A single eigenvalue remains non-zero in the limit 5^0 
(cross). The normalization condition Ai = ~}2 k = 2 corresponds to the rescaling of the width distribution, with (W2) = 1. 



2. Small-window limit (a>3) 

For a > 3, the integral term in eqn (12) is sub-dominant, and the naive expansion of the matrices A pm and A cos in 
powers of 6 becomes correct in the limit 5 — > 0. Using the expansions in eqn (18) and, furthermore, 

C mn = 0(5 4 ), S mn = ^mnir 2 5 2 + 0(5 4 ), I nm = 0(5 3 ), (33) 

we can construct the matrices A per and j4 cos , and check that they have only a single non-zero eigenvalue. For 
illustration, we show in Fig. 7 the spectrum of these two matrices for a = 7/2 at finite S (eigenvalues are normalized 
so that J^k Afe = 2). 

From eqn (23) we then obtain 

«r»(^o) = «n^o) = 

The associated characteristic function reads: 



(21 - 2)!! + 0(5 a - 3 ) fora>3 
(21- 2)!! + 0(1/ log 6) fora = 3 



/(«) = (l-2iAs)-', 
and its inverse Fourier transform gives the universal distribution V(z) for a > 3,: 

exp [—z/2] 



/2ttz 



By comparing the universal distribution for a > 3 with the ones with a < 3 (see Fig. 4), one notices that the additional 
eigenvalues (and branch points contributing to the characteristic function) allow the distribution V(z) to vanish at 
z = 0, thus producing a local maximum in this function. 



IV. CONCLUSIONS 



In conclusion, we have studied the geometric properties of functions with a particularly simple expression in Fourier 
space: independent Gaussian random variables. We have restricted ourselves to one-dimensional self-affine functions 
(characterized by a single length scale), but our analysis evidently carries over to functions with more than one 
length scale, and to higher dimensions. In real space, the geometrical properties are non-trivial, and the boundary 
conditions play an important role. This comes about because the real-space action, for all non-even a, contains 
non-local operators and, for all a ^ 2, is non-Mar kovian. 

We have provided a simple and compact framework for studying the boundary effects for general non-Markovian 
Gaussian processes by relating the characteristic function of the width distribution to the spectrum of boundary- 
dependent infinite matrices, which essentially encode the overlap of the basis functions of the Fourier series. The choice 



12 



of Fourier basis (resulting in the full periodic, cosine and sine series) determines the different boundary conditions. 
We have carried out a complete analysis of the spectrum of these matrices for different values of a for the case of the 
function on the entire interval and also for the function restricted to a window. The associated width distributions 
could all be determined by solving for the eigenvalues of a matrix, and by performing a straightforward inverse Fourier 
transform. 

We have shown that the non-Markovian action propagates the effects of boundary conditions over the entire interval. 
However, in the small-window limit, the width distribution becomes universal (independent of boundary conditions). 
For a < 3, we showed how to compute the cumulants of the width distribution in this limit, avoiding problems related 
to the non-uniform convergence of the Fourier series. For a > 3, the problem of finding the universal width distribution 
drastically simplifies, and we were able to write it down explicitly. Finally, we have obtained the logarithmic corrections 
in the case of odd integer a, in particular for a = 3. 

We hope that our work will be useful for the analysis of experimental data (which usually correspond to our window 
boundary conditions, often in the regime 5 —> 0, which we found to be universal). In many experiments, the roughness 
exponent cannot be extracted reliably by extrapolation, and the width distribution may provide crucial additional 
information. 

Computer programs that compute the width distributions for any value of a, both at finite S and in the limit 5^0, 
are available [27]. 
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APPENDIX A: EXPANSION FORMULA 



In this appendix we derive the expansion formula of eqn (12), first for non-integer, then for integer values of a. 



1. Expansion formula for non-integer a 

We consider the sum 

" f(nS) 



n" 

n—1 

where f(z) is assumed to be a general analytic function. In order to obtain the expansion in powers of 5 of the above 
sum, we expand f(z): 

~ f(n5) " " /(">(()) 

Z^i n a 2^1 2^1 m \ n a-m 

7i—l n—1 m—0 

For a — m < 1, the sum over n converges, and can be directly evaluated using the Riemann Zeta function ((a — m) = 

" f {nS ) ^ /(")(Q)C(m - a) ^ fJnS) 

n—l m—0 n—1 

where [a] is the integer part of a, and f{z) = Y^M f^ m \0)z m /m\. The second term on the right hand side of eqn (Al) 
can be expressed as a contour integral in the complex plane: 

" f(nS) _ f ^ f(z 

71=1 



E n6 o / , J(z fix \ . 
_ V— i d ^ cot (r)< (A2) 



where the contour C encircles the poles z n — nS, n = 1, 2 . . . of the function G(z) = [f(z)/z a ] cot (nz/6). The contour 
C is transformed into the contour C e , as shown in Fig. 8, avoiding the branch point at the origin. By considering the 
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Re z 



Re z 



FIG. 8: Contours for the integral in eqn (A2). The indicated deformation of the contour (C — ► C" — > C e ) transforms eqn (A2) 
into eqn (A3). 



integration over the contour C e , and by performing the change of variables z — ► —\(z/8 — e), we have: 

/[J(is + e)] 



> — — - = -^/ dz ^ coth tt(z - ie ) . 

^ n° 2,/_ 00 iz + e« L V n 



(A3) 



To extract the (5-series of the above integral, it is tempting to integrate term by term the Taylor series of the function 
f(z). We hrst split the above integral into the following terms: 

E 4? = ~\ r dz + £ ^ iz + ^ t COth W z - ^ -n-U^dz f[5(iz + e)](iz + e)-« 

n=l n L Jo L Jo 



dz f[5{iz + e)](iz + e)- a [coth [tt(z - ie)] + 1] + - / dz f[S{iz + e)](iz + e)~ a . 



This separation is valid if all the above integrals are well-defined, a condition respected by the functions f(z) considered 
in this paper (for S < 1). Using the following representation of the Ricmann Zeta function, valid for Re(f3) < 1: 

r o 



C(/3) = lim ^~ dz (iz + e)- /3 (coth [tt(z - ie)] - 1) + J dz (iz + e)-"(cath [tt(z - ie)] + 1) ^ , 



we obtain 



= E Z ^r 1 ^C(a-m) + -^ rdzf{z)(iz + e)-° + U° dz /»(* + e)~ 

„_1 71 r i JO Z J-oo 



(A4) 



Finally, by considering the contours shown in Fig. 9, we verify that 



lim — 

e^O 2 



dz /(<5(iz + e))(iz + e)~ a + / dz f(S(iz + e))(iz + e)~ a = S a ~ 



dt 



Collecting these results proves eqn (12). 



2. Expansion formula for integer a 

The series in eqn (12) is not defined for integer a = m+1, because of the simple pole of the Riemann Zeta function 
((z) at z = 1. However, this divergence is compensated by an ultraviolet divergence in the integral. We analyze this 
situation by considering a = [a\ + p, and by taking the limit p — > 0. We write: 

In order to isolate the ultraviolet divergence, we split the integral into two infrared-divergent terms, and eqn (A5) 
becomes 



J, \ f £j m! I J, (t>-l)»'+' (a -1)1 
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FIG. 9: Contours of the integrals in eqn (A4). 



where e is an infrared cut-off. Using the expansion ((1 + p) = 1/p + 7 + O(p) , where 7 is the Eulcr constant, we can 
take the limit p — > 0: 



lim 



1 



= [7-log((5)-log(e)]. 



From the above we get the following expansion formula for integer values of a: 



(a-l)\ 

where the constant is expressed in term of the following limit: 

'/(*) ^ 2 / (m) (0)i 



C(a — m) 



n=l ^ ' m^a — 1 



const 



V e \ m=0 



m! 



m-a \ r(a-l)(n\ \ 

+ L , U (log(e)+7) ■ 



(a-1)! 



(A6) 



APPENDIX B: GENERAL CUMULANTS 



For completeness, we derive in this appendix the trace formula for the cumulants m of the rescaled width distribution 
\I/(s) for all I from their generating function 



#(s) = log (exp (sio 2 )) = jj-s 



1=1 



We first consider the full periodic series. To compute the average (exp (sic 2)) we have to integrate over the Gaussians 
{a n ,b n }. To make all the Gaussian integrals well-defined, we introduce a cut-off N in the momentum space, such 
that n — 1, . . . , N. The average (exp (SW2)) takes the form: 



(exp(swj )) 



(ll„' I-oc dx "' ) CX P ["I En'ro' M n , ml {s)x n -X v 
(lln> I-oc dx «' ) CX P ["I E„W ^'(Ojln^ 



DetM(s) 



DetM(O) 



where the indices n' and m' run from to 2 AT. In the above expression, we set x n > — a' n , for n' = 0,. . . ,N and 
x n ' — b n /-N for n' = N + 1, . . . , 2 A 7 ". The matrix M(s) is a block of four N x N matrices defined as 



M(s) 











2C n 



2S n 



(for the coefficients C nm , I nm , and S nm , see eqn (9)). We arrive at the following expression for the >F(s): 
*(*) 



■^logDetM(s) + ilogDetM(O) 



iTr[lo g M( S )] + iTr[lo g M(0)]. 
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Expanding the logarithm log M(s) one obtains 



oc 



Tr [(A pcr ) 1 ] 



where A pcr is given in eqn (21). Hence the Ith cumulant of the roughness is given by 



.per 



2(w 2 } 1 



Tr [(A pcr ) 1 ] . 



(Bl) 



For the cosine series, one has to integrate over the Gaussian variables {c„} and {s n }. One arrives at an analogous 
expression for the cumulants. 
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